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We present a method for quantum state tomography that enables the efficient estimation, with 
fixed precision, of any of the matrix elements of the density matrix of a state, provided that the 
states from the basis in which the matrix is written can be efficiently prepared in a controlled 
manner. Furthermore, we show how this algorithm is well suited for quantum process tomography, 
enabling to perform selective and efficient quantum process tomography. 

PACS numbers: QD: 03.65.Wj,03.67.-a,03.67.Pp 



INTRODUCTION 

Quantum information processing tasks always involve 
the preparation and manipulation of quantum systems. 
To be able to perform such tasks it is essential to be 
able to characterize both quantum states and quantum 
operations. The protocols for characterizing a quantum 
state are usually refered to as quantum state tomography 
(QST)[TH5]. In general, QST is a hard task since it in- 
volves an exponentially large number of measurements to 
be preformed (exponentially large on the number of sub- 
systems). Not only that, but the type of measurements 
that one needs to perform on the systems are usually not 
easy to perform. 

On the other hand, the characterization of quan- 
tum processes, known as quantum process tomography 
(QPT)[TJ [7H1], is also an exponentially hard task. How- 
ever, there are some quantum algorithms that allow to 
efficiently extract important information about a given 
quantum process [TUlll6| . These algorithms do not re- 
quire performing QST on the final states but measur- 
ing quantities such as survival probabilities (or transi- 
tion probabilities). Other algorithms for QPT, however, 
do depend upon QST. This makes algorithms for QST 
an essential tool not just for state tomography, but for 
process tomography also. 

In this paper we will present new methods for quantum 
state tomography. To be specific, it is useful to describe 
the quantum states we are to perform tomography in the 
following way:. Let rl be the Hilbert space for the system 
in question, and let B = {\ip a ) > a — !>-••! D} be a basis of 
H, where D — drnxH. Then the density matrix p of a 
state can be written in the basis B as 

D 

P= ^2 Uab\lpa) (lpb\ (1) 
a,b=l 

where a ab = (ip a \p\ipb)- 

In what follows, we will present a method for selec- 
tive efficient quantum state tomography (SEQST) that 



allows to estimate any coefficient a a b with resources scall- 
ing polynomially with the number of subsystems. For 
this to be possible, we need that any state from the basis 
B can be efficently prepared in a controlled manner. This 
method, when applied for implementing QPT results in 
a protocol for efficient and selective QPT that is equiv- 
alent to the one presented in p~3j [14] , illustrating one of 
the virtues of such a selective and efficient QST scheme. 

This paper if organized as follows. First we briefly re- 
view existing methods for QPT that rely on performing 
QST in the final states of a certain process. Then we 
present the selective and efficient algorithm for QST and 
show how it provides the right tool for efficient and se- 
lective QPT, as opposed to previous methods for QST. 
Finally, we compare that QPT algorithm to the one pre- 
sented in [T31 [H] , showing how both can be understood 
in a common theoretical frame. 



QUANTUM PROCESS TOMOGRAPHY BASED 
ON STATE TOMOGRAPHY 

The goal of QPT is to identify the temporal evolu- 
tion enforced by a certain quantum process. Any such 
process is mathematically represented by a linear map £ 
transforming initial states into final states. In fact, the 
operation £ is not only linear but also completely posi- 
tive, and acts as 

£{Pin) = Pout (2) 

This operation represents the discrete (input-output) 
evolution of quantum states. We will focus on maps that 
are trace preserving and whose output dimension is the 
same as the input one. To describe the quantum map it 
is convenient to parametrize it in some way. It is simple 
to notice that any linear map can be written in terms of 
a certain matrix, known as the ^-matrix. This is defined 
with respect to a certain basis of the space of operators. 
In fact, if we choose a basis {E m , m — 0, D 2 — l}, the 
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^-matrix representation for £ is determined by the equa- 
tion: 



£ 



£{P) =^2XmnE m pEl. 



(3) 



This description is completely general since the above 
expressions can be written for any linear channel. Prop- 
erties of the channel arc in one to one correspondence 
with properties of the ^-matrix. In fact, the channel 
preserves the hermiticity if and only if the ^-matrix is 
hermitian. Also, the channel preserves trace, if and only 
if the condition J2mn Xmre^i^m = J is satisfied. Fian- 
lly, the channel is completely positive if and only if the 
^-matrix is positive. Thus, the ^-matrix (which, as we 
mentioned above, depends on the choice of basis for the 
space of operators) fully describes the channel. There- 
fore, quantum process tomography is the task of estimat- 
ing the matrix elements of x- To achieve this goal there 
are several methods, some of which involve performing 
quantum state tomography on final states. Let us review 
them now. 



State Tomography 



FIG. 1: Scheme for the Ancilla Assisted Process Tomography 
quantum algorithm. 



This method, besides requiring n ancilliary qubits, has 
the following two troublesome properties: First, it is not 
clear if the information from the x matrix can be directly 
accessed via measurements on the resulting states. Sec- 
ond but related to the previous point, it is not clear how 
to use QST on the final state to efficiently extract par- 
tial and relevant information on the channel. These two 
issues will be solved in what follows. 



Direct Characterization of Quantum Dynamics 



Ancilla Assisted Quantum Process Tomography 

The Ancilla Assisted Quantum Process Tomography 
(AAPT) E] uses n ancilliary qubits and allows to ex- 
tract all the information about the channel. However, 
as presented in [7J|5], it only allows to obtain full infor- 
mation about the process, being unable to only useful 
partial information about the channel. Thus it is inher- 
ently inefficient as full QPT always is. 

This is the first of two methods that we will men- 
tion here that exploit the state-channel duality given 
by Choi-Jamiolkowski's isomorphism. Such isomorphism 
establishes a one to one relationship between linear op- 
erators from % ®T~L to H <£> H and completely positive 
superoperators acting on the space of operators from % to 
H. The Choi-Jamiolkowski's isomorphisms establishes a 
correspondence between states and channels in the fol- 
lowing way: 

p s =£®n\i)(i\). (4) 

where \I) = ^2j\ii)/y/D is the maximally entangled 
state. After the application of the channel to one of the 
parts, one can perform state tomography to the state ps, 
obtaining full information about the process £. Figure [l] 
ilustrates this algorithm. 

One of the strenghts of AAPT is that the initial state 
can be another state and not necessarily the maximally 
entangled state. As the number of independent parame- 
ters defining the initial state (the Shmidt number) is D 2 , 
such state can always encode the necessary information 
to define the quantum channel. 



The Direct Characterization of Quantum Dynamics 
(DQCD)[H1 [5] is a quantum algorithm similar to that 
of AAPT in many aspects. It also resorts to n ancil- 
liary qubits, and relies on the Choi-Jamiolkowski iso- 
morphism. Contrary to the AAPT, on the DCQD the au- 
thors explicitely showed a method to efficiently and selec- 
tively measure the diagonal coefficients of the ^-matrix; 
however, off-diagonal elements still require to invert en 
exponentially large matrix. This makes the method inef- 
ficient. 

To describe the method let us consider the operator 
basis consisting of the n-fold tensor product Pauli oper- 
ators acting on each qubit. We denote these operators 
as \E m ,m — 0, ...,D 2 — l}. In that basis, the channel 
description is given by 

£(p) = ^2 XmnEmpEt ( 5 ) 
mn 

In order to perform diagonal tomography, the algo- 
rithm proceeds as follows. First, as the AAPT, we have 
to generate the state that is isomorphic to the channel 

P£ = h^ XmnE m \i)(j\El ® (6) 

ijmn 

Now, the probability of measuring the state J2i 
on the output is given by 

i J2(kk\ P£ \U) = ^ E X«™ Tr Tr ( E n) = Xoo. 

kl mn 

(7) 

Thus, we see that the survival probability of the input 
state directly gives the coefficient xoo- It is easy to show 
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that the very same method can be used to evaluate any 
diagonal coefficient of the x matrix. In fact, the prob- 



ability of obtaining the final state 
nothing but Xkk- As the set 



J2i E k ® l\H) is 



TZ= {E k (Z>I\I),k = 0,...,D 2 -1} 



(8) 



forms an orthonormal basis of % ® H, a measurement 
in that basis will suffice for diagonal tomography of the 
%-matrix. 

The main problem arises when off-diagonal tomogra- 
phy is taken into account. The solution presented in [8l|9] 
is to use a state other than a maximally entangled one. 
However, it can be shown that in the most general case, 
this procedure requires inverting an exponentally large 
matrix. Again, this makes the method inefficient. 

Below, we will introduce a method for QST that pro- 
vides not only a convenient tool for QST, but also would 
provide the necessary ingredient missing in AAPT and 
DCQD to obtain an efficient and selective QPT. 
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FIG. 2: Quantum circuit for Selective and Efficient Quantum 
State Tomography. 



basis B = {\ip a ), a — 1, D} from the state \ipo). That 

IS, V a \i>o) = \lpa)- 

It is straightforward to verify that by measuring the 
average value of the operator |-0o) (V'o I <8 cr x (that is the 
average value of the operator a x conditioned on the de- 
tection of the state |?/>o) on the main system) one obtains 
the real part of Xab- Moreover, replacing a x by a y , the 
same method provides the imaginary part of the same 
matrix element. Thus, 



SELECTIVE AND EFFICIENT QUANTUM 
STATE TOMOGRAPHY 



Tr(p F \*p )(ip \®a x ) = SRxafc (11) 
Tr(p F \ipo){ipo\®<T v ) = 3 X afc (12) 



The standard method for QST was clearly described in 
P]. This method resorts to the description of the state 
in the Pauli operator basis as 



(9) 



where Ei are the n-fold Pauli operator basis for an n 
qubit system. It is straightforward to perform tomogra- 
phy in this basis by just measuring the expectation value 
of every Ei. Although this method is indeed selective 
and efficient, it is not well suited for selective and effi- 
cient QPT. 

The method we are about to introduce, the Selective 
and Efficient Quantum State Tomography (SEQST), is 
also efficient and selective but, as opposed to the stan- 
dard method, it is selective in any basis of the corre- 
sponding Hilbert space. That is, given the state written 
in the form 



D 



(10) 



a,b=l 



and provided we know how to prepare the states from 
the corresponding basis in a controlled manner, we will 
be able to selectively measure any given coefficient a a f,. 
Such measurement will be efficient, meaning that given 
a precision, the number of required single click measure- 
ments doesn't scale with the size of the system. 

Consider the circuit shown in Fig. [2] where the oper- 
ators V a are the ones that prepare the states from the 



where pp is the state prior to the measurement. 

To discuss the efficiency of the method we should an- 
alyze the resources needed by this algorithm. First of 
all, the efficiency of the method is limited by that of the 
implementation of the controlled-Vj and controlled-Vj 
operators. In fact, if the implementation of such oper- 
ations require O (f (n)) operations, then the full circuit 
will also require O (/ (n)). 

The other point to determine the efficiency of the 
method is to determine the number of experimental runs 
required to obtain the desired result to a given precision 
e with a probability of success p. To answer that ques- 
tion we just need to consider that each experimental run 
gives one of three results (+1 corresponding to |^o) on 
the system and |0) on the ancilla, — 1 corresponding to 
IV'o) on the- system and |1) on the ancilla, and corre- 
sponding to another result on the system). The % ma- 
trix element is estimated by computing the average value 
these results after performing a certain number of repeti- 
tions M. Each of the results are detected at random with 
their corresponding probabilities. Therefore, one can use 
a Chernoff bound to show that to obtain the correct re- 
sult with uncertainty e and a probability p of success, the 
number of experimental runs M must be such that 



2 In 



M > 



(13) 



which does not depend on n or D. This implies that the 
method is efficient. 
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Application to Quantum Process Tomography 

In this section we will show how the SEQST algorithm 
is the right tool to perform selective and efficient QPT 
when combined with the AAPT method reviewed in Sec. 

D 

In order to proceed, we need to find out the way in 
which the quantum state isomorphic to the channel £ 
depends on the ^-matrix of such channel. Thus, we will 
show that the ^-matrix of the channel is nothing but the 
matrix element of the quantum state in a particular basis. 
Therefore, by performing quantum state tomography in 
that basis we directly provide the information about the 
%-matrix of the channel. To show this, we use equation 
Q and replace the expansion of £ in the Pauli basis: 



PS =^2 XmnE m ® I|i) (I\E% <g> I 



(14) 



Indeed, this shows that x is the matrix representation of 
in the basis 1Z shown in equation (JsJ) . Therefore, to 
selectively measure a single \ coefficient one only needs 
to perform selective tomography in the basis 1Z. 

Using the results already presented in the previous sec- 
tion, we see that to do this we should implement the cir- 
cuit described in Figure [3] Here, the application of the 
channel £ to one of the pieces of the maximally entan- 
gled state can be regarded as the preparation of the state 
isomoprhic to the channel. In turn, the rest of the circuit 
is nothing but the SEQST algorithm described above. 
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FIG. 3: Application of SEQST to QPT. 



It is important to point out that, since the measure- 
ment is direct, the analysis of the resources required to 
implement the method (presented in the previous sec- 
tion) directly applies to this case. The only extra re- 
sources needed in this case are involved in the prepara- 
tion and measurement of the maximally entangled state 
which require O (n) aditional single and two qubit gates. 



COMPARISON WITH OTHER SELECTIVE AND 
EFFICIENT QPT SCHEMES 

Another quantum algorithm for selective and efficient 
quantum process tomography is the one known as SE- 
QPT, precisely for Selective and Efficient Quantum Pro- 
cess Tomography[T31[T3]. The main idea there is to follow 



the proceedure described by the circuit shown in Figure 
|4j and to estimate the average answer averaging over the 
entire Hilbert space of the system using the Haar mea- 
sure for that purpose. As it is shown in [T31 Q3], that 
average cn be directly related to the matrix element \ab 
as 



(£r x ®|V)<## = 
(a y ® \ip){tp\)d^ = 



£>RC (Xab) + $ab 

D+l 
glm (Xab) 
D + l 



(15) 
(16) 



Moreover, it can be shown that the average over the en- 
tire Hilbert space can be efficiently estimated by ran- 
domly sampling over a special (and finite) set of states 
which is known as a 2-designs. For these reasons, the 
method SEQPT is not only selective but also efficient. 
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FIG. 4: Circuit for the Selective and Efficient Quantum Pro- 
cess Tomography algorithm. Depending on the measurement 
of (j x or <Jy, the real or imaginary part of Xafo will be obtained. 

As we mentioned above, in the SEQPT scheme, the 
average is estimated by randomly sampling states. In 
the scheme we proposed above, the average is obtained 
automatically by the quantum correlation between both 
parts of the maximally entangled state, without the need 
to resort to randomly preparing and detecting the special 
states of the 2-design. 



SUMMARY 

In this paper we presented a novel quantum algorithm 
to perform selective and efficient quantum state tomog- 
raphy in any Hilbert space basis, given that the states 
from that basis can be efficiently prepared in a controlled 
manner. 

We then showed that, when properly combined with 
the Ancilla Assisted Process Tomography, it yields a pro- 
tocol for QPT that is both selective and efficient. Finally, 
we showed that this protocol shares some properties with 
SEQPT, a method presented in [13 EJ. The main dif- 
ference is that the use of ancillas is a way to avoid the 
preparation of the special states of the 2-design and sam- 
pling on them. 

This work was partially supported by funds from CON- 
ICET, ANPCyT, UBACyT, EU Q-Essence and Spanish 
FIS2010-14830 projects. A.B. is funded by an ERC Start- 
ing Grant PERCENT. 



5 



[1] M. A. Nielsen and I. L. Chuang, Quantum Computation 
and Quantum Information (Cambridge University Press, 
2000). 

[2] M. Paris and J. Rehacek, Quantum state estimation, Lec- 
ture notes in physics (Springer, 2004) 

[3] D. Gross, Y.-K. Liu, S. T. Flammia, S. Becker, and 
J. Eisert, Phys. Rev. Lett. 105, 150401 (2010) 

[4] M. Cramer, M. Plenio, S. Flammia, R. Somma, D. Gross, 
S. Bartlett, O. Landon-Cardinal, D. Poulin, and Y.-K. 
Liu, Nature Co mm. 1, 149 (2010) , [http : //arxiv . org/ 
| abs/1101. 43661 1101.4366 

[5] C. Miquel, J. P. Paz, M. Saraceno, E. Knill, R. Laflamme, 
and C. Negrevergne, Nature 418, 59 (2002) 

[6] A. Bendersky, J. P. Paz, and M. T. Cunha, Phys. Rev. 
Lett. 103, 040404 (2009) 

[7] J. B. Altepeter, D. Branning, E. Jeffrey, T. C. Wei, P. G. 
Kwiat, R. T. Thew, J. L. O'Brien, M. A. Nielsen, and 
A. G. White, Phys. Rev. Lett. 90, 193601 (2003) 



[8] M. Mohseni, A. T. Rezakhani, and D. A. Lidar, Phys. 

Rev. A 77, 032322 (2008) 
[9] M. Mohseni and D. A. Lidar, Phys. Rev. Lett. 97, 170501 

(2006) 

[10] J. Emerson, M. Silva, O. Moussa, C. Ryan, M. Laforest, 
J. Baugh, D. G. Cory, and R. Laflamme, Science 317, 
1893 (2007) 

[11] C. C. Lopez, B. Levi, and D. G. Cory, Phys. Rev. A 79, 
042328 (2009) 

[12] B. Levi, C. C. Lopez, J. Emerson, and D. G. Cory, Phys. 

Rev. A 75, 022314 (2007) 
[13] A. Bendersky, F. Pastawski, and J. P. Paz, Phys. Rev. 

Lett. 100, 190403 (2008) 
[14] A. Bendersky, F. Pastawski, and J. P. Paz, Phys. Rev. 

A 80, 032116 (2009) 
[15] C. T. Schmiegelow, A. Bendersky, M. A. Larotonda, and 

J. P. Paz, Phys. Rev. Lett. 107, 100502 (2011) 
[16] C. C. Lopez, A. Bendersky, J. P. Paz, and D. G. Cory, 

Phys. Rev. A 81, 062113 (2010) 



